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We show that the width of the wave-packet of a class of 
generalized nonlinear Schrodinger equations (NLSE) trapped 
in an arbitrary time-dependent harmonic well in any dimen- 
sions is universally determined by the same Hill's equation. 
This class of generalized NLSE is characterized by a dynam- 
ical 0(2, 1) symmetry in absence of the trap. As an appli- 
cation, we study the dynamical instabilities of the rotating 
as well as non-rotating Bose-Einstein condensates in one and 
two dimensions. We also show exact extended parametric res- 
onance in a non-relativistic Chern-Simons theory producing 
a gauged NLSE. 

PACS numbers: 05.45.Yv, 03.75.Fi, 11.15.-q, 03.65.Ge 



The nonlinear Schrodinger equation (NLSE) appears 
in many branches of present-day physics and mathemat- 
ics The optical soliton a solution of the NLSE, 
has even been observed experimentally 0. The one di- 
mensional NLSE is exactly solvable. There exists several 
other generalized NLSE in one dimension which are also 
exactly solvable. However, a very little exact and analyt- 
ical results are known for higher dimensional generaliza- 
tions of these models, although they are very much rele- 
vant in many branches of modern science. The purpose of 
this letter is to present an exact, analytical description of 
the dynamics of the width of the wave-packet for a class 
of generalized NLSE in arbitrary dimensions trapped in a 
time-dependent harmonic well. This class of generalized 
NLSE is characterized by a dynamical 0(2, 1) symmetry 
in absence of the trap. 

Consider the following Lagrangian in arbitrary d + 1 
dimensions. 



1+7*' 



1 

2m 



| VV> \-gV(il>,tp*,T). 



(1) 



The coupling constant g has the inverse-mass dimension 
in the natural units with c = h = 1 . The real potential V 
does not depend on any dimensional coupling constant. 
This allows to have a scale-invariant theory. We demand 
the invariance of the action A = J dTd d rC under the 
following time-dependent transformations 



r -> r h = f(t) ^ r , T _> t = t(r), f(t) 



dr(t) 
dt ' 



ip(r,T) -> iph(t,r h ) = T^exp ( -im—r h ) t/>(r,r), (2) 



with the scale-factor r given by. 

at + (3 



Particular choices of r(t) = t + f3, a 2 t, and 
spond to time translation, dilation and special conformal 
transformation. The generators of these transformations 
( as given below ) close under an 0(2, 1) algebra. Al- 
though V is restricted to have specific forms due to the 
requirement of the symmetry, one can still make infinitely 
many choices of it. ft might be noted here, apart from its 
dependence on the condensate ip, the potential can also 
be explicitly dependent on the space-coordinates. We 
keep V arbitrary, but, consistent with the 0(2, 1) sym- 
metry in ([!]), unless mentioned otherwise. 

Let us now introduce two moments I\ and I2 in terms 
of the density p and the current j as, 

p(r,r) =tp*ip, j(r,r) = --L(^*V^-VV^*), 
h(r) = ^jd d v r 2 p, h{r) = | J d d v r • j. (4) 

We are dealing with a conservative system and the to- 
tal number of particles N(t) = J d d rp is a constant of 
motion. The moment I2 is related to the speed of the 
growth of the condensate. The moment I\ describes the 
square of the width of the wave-packet H. This quan- 
tity plays the central role in the analysis of the collapse 
of the condensates of the NLSE with or without time- 
independent harmonic trap [fz| — [lO^I - ft is also used in the 
context of the Bose-Einstein condensation(BEC) |llf to 
study the low energy excitations and in optics |12j] to de- 
termine the beam parameter evolution. The dynamics of 
Ii , when the system (|]) is immersed in an external time- 
dependent harmonic trap, is the central subject of the 
investigation of this letter. We show that the dynamics 
of 1\ is universally determined by the same solvable Hill's 
equation, independent of the space-dimensionality, inte- 
grability and nature ( short-range, long-range, local, non- 
local, linear, non-linear ) of the interaction. The univer- 
sality in the description of the dynamics of the width for 
this class of theory has been observed partially through 
a time-dependent variational analysis in |L3| . We present 
here exact, analytical and complete treatment. 

The system ([j]) has a dynamical (9(2, 1) symmetry. 
The generators, the Hamiltonian H, the dilatation gen- 
erator D and the conformal generator K are, 



H = d d r 



-L\V^\ 2 +gV(^,r,r) 
Ira 



r(i) 



7* + S ' 



aS — /?7 = 1. 



(3) 



D = tH-I 2 , 

K = —t 2 H + 2tD + I\ . 



(5) 
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These generators are constant in time and lead to the 
following equations, 



dH dL dl 2 

— = 0, -^ = 2/2, -r l = H. 
dr dr dr 



(0) 



For time independent solutions, both I\ and I2 do not 
depend on r. As a consequence, the static solutions of 
a system with 0(2, 1) symmetry carry zero energy [ fl4"|| . 
We also note that I2 = D = and K — 1\ for static 
solutions of Defining the width of the wave-packet, 
X = \/Ti, it is easy to find a decoupled equation for X 
fromfl), 



d 2 X _ Q_ 
~d^~X^' 



Q = hH 



Io > 0, 



dQ 
dr 



0. 



(7) 



The constant of motion Q is the Casimir operator of the 
0(2, 1) symmetry. Eq. (|t]) can be interpreted as the 
equation of motion of a particle moving in an inverse- 
square potential. Interestingly enough, this system also 
has a dynamical 0(2, 1) symmetry. This reduced sys- 
tem of a particle in an inverse-square potential is a well- 
studied problem and the solution is given by [Q , 



X 2 = {a + b T f + %T\ 



(8) 



where a and b are the integration constants. 

Consider the time-dependent transformations in (|^) 
with arbitrary scale- factor r(t). This is no more sym- 
metry transformations of (Q) for general r. The action 
A is transformed into a new one, Ah = j dtd d rhCh, con- 
taining a time-dependent harmonic trap. The new La- 
grangian Ch now reads, 



C h = iip* h dttph 



2m 



gV(ip h ,ipl,r h ) - -mu(t)rl\ ip h | 2 . 



(9) 



The time-dependent frequency u>(t) of the harmonic trap 
is determined by, 



b + uj(t)b = {), b(t) = r 



(10) 



Once the solution of the equation of motion of ([[]) is 
known, the same can be obtained for (H) by using the 
transformation (0) and the equation ( J10[ ) or the vice 
versa. Equation (nu) describes the motion of a particle in 
a time-dependent harmonic trap. For r(t) — —tan(uiot) , 
it gives u> = ujq. For the special choice of the fre- 
quency u) obviously vanishes. The general solution of jlO| ) 
for the physically relevant periodic u(t) is well-known and 
will be discussed below. 

The dynamical 0(2, 1) symmetry of C is not present 
for Cu. We replace (r,ip,r) by (rh,iph,t) in the defini- 
tion of I\ , I2 and H and denote the resulting expressions 
in terms of the 'curly form' of the associated variables. 



Under the transformation 
following form, 



), the equations (a) have the 



i 1 (t) = 2i 2 {t), 
i 2 {t)=H(t)-cj(t)i 1 {t) 7 

H(t) = -2uj(t)I 2 (t). 



(11) 



Defining a new variable X(t) = ^/Ti(t), it is easy to find 
a decoupled equation for X ', 



X+co(t)X=^, 
Q=l 1 H-li>0, Q = 0. 



(12) 



We have the surprising result that the dynamics of the 
width of the wave-packet of the system (|^) is universally 
determined by the equation (^). This result is inde- 
pendent of the integrability of the model. We also have 
the freedom of choosing a large class of V as long as the 
dynamical 0(2, 1) symmetry in absence of the harmonic 
trap is maintained. The knowledge of the time-evolution 
of X allows us to determine the time-evolution of H, 



n = x 2 



Q_ 
x 2 ' 



(13) 



For the time-independent trap, ui(t) = loq, the Hamilto- 
nian Hh = H + uqT\ corresponding to the Lagrangian 
Ch is a constant of motion. The Hamiltonian Hh is re- 
lated to the generator of the compact 5*0(2) rotation of 
£0(2, IV 

Eq. ( |l2| ) can be interpreted as that of a particle moving 
in a time-dependent harmonic trap and a inverse-square 
potential. Due to the underlying 0(2, 1) symmetry in ab- 
sence of the trap, this equation can be obtained directly 
from (Q) through the use of the transformations (Q) . The 
dynamics of the width X can thus be constructed exactly 
from d|), 



X(t) = b(t)X( T (t)), 



(14) 



with the knowledge of the scale-factors r(t) and b(t) from 
( pTj| ) for a particular choice of u>(t). We provide below a 
familiar form of solution of Eq. ( |l2| ) , since it appears in 
many branches of physics including the cylindrically sym- 
metric two-dimensional NLSE |TJJ] . The general solution 
of (|T^) is given by, 



X 2 (t)=u 2 {t) + ^v 2 {t), W(t) 



(15) 



where u(t) and v(t) are two independent solutions of the 
equation, 



X + Uj(t)x = 0, 



(16) 



satisfying u(to) = X(to),ii(t ) — X(to),v{to) = and 
v(t ) ^ 0. For periodic ui(t) with the period T, the above 
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equation is known as the Hill's equation and is a text- 
book material |51|. We just mention here the general sta- 
bility criteria in terms of the quantity 5 — \ u(T) + v(T) \ 
with the normalization X(0) = 0, X(0) = 1 and v(0) = 1. 
The solution is stable for 5 < 2, while it is unstable for 
5 > 2. We remark that the same stability criteria is valid 
for Eq. M). 



The central result of this letter is contained in (12) 



Although our main concern in this letter is on NLSE, we 
remark that the same result is true for a class of linear 
Schrodinger equations with Calogero-type inverse-square 
interaction in arbitrary dimension |16|. An example of 



arbitrary d dimensional non-linear potential consistent 
with 0(2, 1) symmetry is given by, 



V(^,r,r) = J d d r>*(r')C/(r-r')^(r')| ^(r) \\ (17) 

with U(r) having the following scaling property. For r — > 
er, U(r) — > U(er) — e~ 2 [/(r). We now discuss a few 
specific examples of NLSE with different choices of V 
which are relevant in the contemporary literature. 

BEC ind=l,g>0,V =\ipf : The Gross-Pitaevskii 
equation (GPE) describing the repulsive Bose-Einstein 
condensates trapped in a time-dependent harmonic trap 
in one dimension can be obtained from the Lagrangian 
® WA - Exact soliton solutions of the GPE equation 
have been obtained in absence of the trap 17 lq]. Only 



approximate or numerical results are known, when the 
time-independent harmonic trap is included fllT]]. 

Following our analysis, the exact solutions of (|9j) can be 
obtained from those of (|]) by simply using the transfor- 
mation (0) and the equation (|io| ) determining the time- 
dependent scale-factor r for a particular choice of the 
Lu(t). We consider the case of time-independent trap with 
oj(t) = ujq and choose g = Define the following di- 
mensionless variables, 



Xfl = TTlp Q X h , t 



7T^ 



t, tph 



mU} ° MS) 



where ipQ, the asymptotic value of the density, is related 
to the chemical potential fi by, ip^ = V 2ot/j . The exact 
solution for iph is, 



iph 



1 



7=~J\ exp 

cosyujQt) 



tan[u>( S t)Xf l — —taniuo^t) 

2 ujq 



("0*) 



- 1 



cosh[- 



Jv_ 



(19) 



In the limit luq — > 0, the solution for the system without 
the trap is recovered [[l7| . Without loss of any generality, 
we are choosing (3 = in Eq. (12) of |l7[ . 

BEC in d = 2, V = \ if) f: We get the GPE describing 
two dimensional BEC in a time-dependent trap from the 
Lagrangian (||) . No exact solution of this GPE with or 



without the trap is known. However, this is a well-studied 
system and many of the dynamical properties are already 
known |T^,^,|lO|. We concentrate here on the rotating 
BEC and present some new interesting results. Consider 
a further time-dependent rotation in (|^) 0] , 

r> = (™f® Sm {SW (20) 
■smf(t) cosf{t) I 



t — > t = t, Th 



This transforms Ah to Ah = j dtd d VhCh containing an 
additional term proportional to the z component of the 
angular momentum with the coefficient given by a time- 
dependent frequency. In particular, the new Ch is given 
by, 

1 2 

C h = itph d i^h - 7T- 1 V/jVft. I - gV(ip h , ip* h , r h ) 
Zm 

1 



~rmj(t)rl\ ip h \ - f ip* h L z ->p h , 

L z = -l\ XhTTZ Vh-^Z- 

Oljh ox h 



(21) 



where Xh and yh are the components of the two dimen- 
sional vector ih- This is the Lagrangian for rotating 
BEC in an external time-dependent isotropic trap in two 
dimensions Interestingly, once the solution of the 

equation of motion of (|l|) is known, the same can be ob- 
tained for (|2l|), using the equations (|2|), JlC| ) and (2C), 
or vice versa. Moreover, under the transformation (2C), 
the set of equations in (|ll]) remains the same in terms 
of the new variables t and ?/,. Thus, the dynamics of 
the width X(t) of (^l|) is again universally determined 
by the equation (|l^). The introduction of the last term 
in ( ^l| ) does not change the dynamical properties of the 
width. There may be dynamic instabilities solely due to 
the rapid fluctuations in the phase of the condensate dur- 
ing the evolution in time. However, it is obvious from the 
definition of I\ that such instabilities do not show up in 
the evolution of the width. 

Gauged NLSE: We now show that the dynamics of the 
width remains unchanged even if the non-trivial gauge- 
fields are introduced in ([j]) maintaining the 0(2, 1) sym- 
metry. Consider a Lagrangian in 2 + 1 dimensions with 
the gauge-fields (Aq, A) and the matter-field if>, 



L g = iiP*{d T - iA )^ - _L| (V - zA)<A | 2 
Zm 



flVir.r . v] +± e ^ Fllv A x , 



(22) 



where the last term is the Chern-Simons term. The mo- 
ment I\ for this non-relativistic Chern-Simons(CS) the- 
ory can be interpreted as the width of the soliton or al- 
ternatively as the quadrupole moment. For V — ^| if) | 4 , 
this is the Jackiw-Pi model describing gauged NLSE [g2| . 
This is relevant in theories with anyons and in the quan- 
tum Hall effect |f22|-^4|j . The Jackiw-Pi model is exactly 
solvable at the self-dual point, g = ^j^y- Our result is 



3 



valid for arbitrary V maintaining 0(2, 1) symmetry. For 
the particular case of Jackiw-Pi model, the importance 
of our result lies at all non-self-dual points, where the 
model is not integrable. The Hamiltonian is given by, 



H = d 2 r 



2m 



(23) 



The CS term being a topological term do not contribute 
to the Hamiltonian. The generators D and K have the 
same expressions as in (|^), with the partial derivative in 
the expression of the current j in the definition of the mo- 
ment I2 replaced by the respective covariant derivative. 
The treatment is now identical to the case without the 
gauge-fields. The same transformations (^) with d = 2 
and the gauge-fields transforming accordingly, 



x" = (T,r), A» = (A ,A), = 0,1,2, 



(24) 



introduce a time-dependent harmonic trap 22 24 1 . The 



width of the soliton of this new Lagrangian is again uni- 
versally determined by (12). Interestingly, the introduc- 
tion of the gauge-fields and the nontrivial CS term to the 
usual two dimensional NLSE does not change the dy- 
namics of A". A comment is in order at this point. It is 
known that Eq.(|i~2|) admits parametric resonances. Thus, 
the solitons of the non-relativistic CS theory should ex- 
hibit the same phenomenon. This provides an example of 
exact, extended parametric resonance in a gauge theory 
with the nontrivial CS term. 

In conclusion, we have shown that the width of the 
wave-packet of a class of generalized NLSE is universally 
determined by the same Hill's equation. This class of 
NLSE is characterized by a dynamical 0(2, 1) symmetry 
in absence of the trap. The result is so robust that it is 
independent of, (i) the space dimensionality, (ii) the inte- 
grability of the model, and (iii) short-range, long-range, 
local, non-local, linear or non-linear nature of the many- 
body interaction. This result persists with its full gener- 
ality even when the gauge-fields are introduced maintain- 
ing the dynamical 0(2, 1) symmetry. The later example 
allows us to study an exact parametric resonance in a 
theory with the nontrivial gauge fields. Special cases of 
this class of generalized NLSE are relevant in BEC and 
in non-relativistic Chern-Simons theory. It would be nice 
to see the importance of this class of NLSE in many more 
physical systems. 
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